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Formation Gibbs free energy of point defects and oxygen clusters in uranium dioxide at high- 
pressure high-temperature conditions are calculated from first principles, using the LSDA-I- U ap- 
proach for the electronic structure and the Debye model for the lattice vibrations. The phonon 
contribution on Frenkel pairs is found to be notable, whereas it is negligible for the Schottky defect. 
Hydrostatic compression changes the formation energies drastically, making defect concentrations 
depend more sensitively on pressure. Calculations show that, if no oxygen clusters are considered, 
uranium vacancy becomes predominant in overstoichiometric UO2 with the aid of the contribution 
from lattice vibrations, while compression favors oxygen defects and suppresses uranium vacancy 
greatly. At ambient pressure, however, the experimental observation of predominant oxygen defects 
in this regime can be reproduced only in a form of cuboctahedral clusters, underlining the impor- 
tance of defect clustering in \J02+x. Making use of the point defect model, an equation of state 
for non-stoichiometric oxides is established, which is then applied to describe the shock Hugoniot of 
UO2+1. Furthermore, the oxidization and compression behavior of uranium monoxide, triuranium 
octoxide, uranium trioxide, and a series of defective UO2 at zero Kelvin are investigated. The evo- 
lution of mechanical properties and electronic structures with an increase of the oxidation degree 
are analyzed, revealing the transition of the groundstate of uranium oxides from metallic to Mott 
insulator and then to charge-transfer insulator due to the interplay of strongly correlated effects of 
5/ orbitals and the shift of electrons from uranium to oxygen atoms. 
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I. INTRODUCTION 

Uranium oxides, especially uranium dioxide (UO2), are interesting materials from a fundamental point of view, not 
only because of their complex electronic structure arisen from partially filled 5/ electron shells and as materials of 
great technological importance in nuclear applications, but also because they are typical nonstoichiometric compounds. 
UO2 manifests fluorite structure at all temperatures up to the melting point and pressures up to 40 GPafiiS as well as 
a wide range of off-stoichiometry.— A comprehensive understanding of its thermodynamical, structural, and kinetic 
properties under extreme conditions is very important, especially for nuclear industry since UO2 is the main fuel 
component for pressurized water reactors. Point defects are the key contents for these properties, particularly under 
irradiation as they play a major role in atomic diffusion. They also act as accommodation of the strong stoichiometry 
variations that exist in this material and provide incorporation sites for the fission products. 

There are some theoretical studies on the properties of point defects^"— and their incorporation with fission 
productsi^"— in UO2 at atomic-scale. Early literature used mainly empirical potentials j^iii while recently the 
ones based on density functional theory (DFT) prevail. Since the standard local density approximation (LDA) or 
generalized gradient approximation (GGA) for the exchange-correlation functional fails to reproduce the fact that 
UO2 is a Mott-Hubbard insulator with an energy gap split by the strongly correlated uranium 5/ electrons^^"— 
current state of the art approaches make use of DFT-I- U formalism or the hybrid functional method to describe the 
electronic structure of U02i^'^"— With these methods, the formation energy of point defects in UO2 at zero Kelvin 
were extensively studied,— including both neutral and charged defects. The population of point defects at 
finite temperatures and partial oxygen pressures were also investigated based on these K formation energies. 

After these efforts, understanding of the electronic structure, magnetic ordering, atomic structural geometry, and 
energetics of perfect UO2 and single point defects in it was advanced enormously. Nevertheless a fundamental 
experiment observation that oxygen defects predominate in the hyperstoichiometric regime still can not be reproduced, 
even at a qualitative leveli^^i^^ Only one LSDA-I- U calculation predicted the predominance of oxygen interstitial in this 
regionjS, However, the agreement is quite marginal. On the other hand, if taking oxygen defect clustering into account, 
the predominance of oxygen defects in a form of cuboctahedral cluster can be easily establishedf21i2S indicating the 
complicated nature of this problem. 

Furthermore, almost all of the previous investigations predicted a negative formation energy for single oxygen 
interstitial Oi, except one of the 'DFT+U calculations with a GGA exchange-correlation functional that predicted 
a positive value of 0.1 eVi^^ This is sharply in contrast to other DFT+U and empirical results which usually are 
about — 1 — 3eVi^iiSii^ It does not sound unreasonable considering that the hyperstoichiometric phase of UO2 is 
unstable at OK.^- However, it is worthwhile to note that a negative formation energy of O^ does not conflict with this 
observation too, since oxygen clusters have a lower energy2i21i^i^ and phase separation between the ordered phases 
of fluorite UO2 and 11403^"— could also lead to the same phase diagram. Nevertheless, a formation energy of Oi 
being apparently positive indicates a different oxidization mechanism, in which the lattice vibrations, or equivalently 
the phonons, should play a pivotal role at finite temperature and therefore alter the sign of the formation energy. This 
part of the contribution to the stability of defects in U02±3;, however, was totally absent in almost all of the previous 
investigations. It now becomes imperative to assess the role played by this part of free energy on the energetics of 
defects. 

In addition to reducing the free energy, lattice vibration at high temperature also leads to thermal expansion. The 
resultant dilation of the system volume inevitably changes the equilibrium K formation energy. The influence of 
this effect on defect populations, as well as that due to possible residual stress arisen from irradiation damages, are 
unclear, which is possibly a partial consequence of the fact that the hydrostatic compression behavior and the equation 
of state (EOS) of solid U02±a; drew little attention in the past decadesi^ These mechanical properties now become 
important. It is crucial to know how the fuel components respond to shock waves generated at nearby explosions 
to assess the general safety of nuclear reactors and spent fuel storage facilities. Attempts to address these problems 
require a profound understanding of defects behavior at high-temperature high-pressure conditions, as well as a full 
EOS for non-stoichiometric fuel compounds. 

In this paper, we evaluate the importance of lattice vibration and hydrostatic compression on the energetics and 
population of point defects in UO2 with DFT-I- U calculations. An equation of state model for off-stoichiometric U02±a: 
is then developed. It will benefit the understanding of other compounds since UO2 is a typical non-stoichiometric 
oxide, and an EOS theory developed for it can be easily adapted for other alike materials. We also attempt to answer 
the question of whether oxygen clusters are definitely necessary or not in order to make oxygen defects the major 
component when contributions of phonon and thermal expansion have been taken into account. In addition, in order 
to demonstrate the overall trend of oxidation effects and highlight the crucial role played by U02±2;, the energetics 
and electronic structures of uranium monoxide (UO), triuranium octoxide (UaOs), and uranium trioxide (UO3) are 
investigated as well. The calculation methodology and the theoretical formalism are presented in the next section. 
The results and discussion are given in Sec. HI, which is then followed by a summary of the main conclusions. 
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TABLE I. First-principles results for the energy curves of UO, \J02±x, Q-UsOg, and 5-UO3, where x is the deviation from the 
stoichiometric composition of uranium dioxide and A'^ is the total number of atoms in the simulation cell. D (in eV), ro (in A), 
and Bq (in GPa) are the cohesive energy per atom, the equilibrium lattice parameter of the effective cubic cell, and the zero 
pressure bulk modulus, respectively. 
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II. METHODOLOGY 
A. Ground-state energy calculation 

The energetics of UO, U3O8, UO3, and a series of defective UO2 in fluorite structures (to model off-stoichiometric 
U02±3;) were studied through calculating the total energy curves with the plane-wave method using the density 
functional theory to treat the electronic energy as implemented in the VASP code^S^^^ The electron-ion interac- 
tions were described by projector-augmented wave (PAW) pseudopotentialsi^i^ The 2s^2p^ electrons in oxygen and 
6s^6p^5f^6d}7s^ in uranium were treated in valence space. For UO2, the electronic exchange-correlation energy was 
computed by spin-polarized local density approximation with an effective on-site Coulomb repulsive interaction to 
split the partially filled 5/ bands that localized on uranium atoms (LSDA+ [/)j^"— The LSDA+C/ energy functional 
is given by 

E-LSDA+U = ^'lSDA + ^'Hub — -E^dc, (1) 

where i^LSDA is the LSDA contribution to the energy, E'nub is the electron-electron interactions from the Hubbard 
term, and i?dc is the double-counting correction. The LSDA+ U approximation is thus a correction to the standard 
LSDA energy functional. Enuh and E^c depend on the occupation matrices of the correlated orbitals. We used the 
simplified rotationally invariant approximation for the Hubbard term i?Hub due to Dudarev et al^ As regards the 
double counting correction, we used the fully localized limit, whose expression is 

E^. = ^N{N - 1) - ^ ^ 7V-(iV- - 1). (2) 

(7 

Parameters of the Hubbard term were taken as U = 4.5 eV and J = 0.51 eV, which have been checked carefully 
for fluorite U02'^'^'^>^ The possible metastable states introduced by strong correlation effects were removed via 
a quasiannealing procedure,^i which has the capability of yielding the same results as an alternative but powerful 
approach of monitoring the occupation matrices . ^^''^'^ Previous investigation implied that the value of the Hubbard 
parameters depends on the local coordination environment,^ and thus it should not be the same for different oxides. 
Since there are no appropriate values of U and J available for UO, U3O8, and UO3, we used the pure GGA with 
the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional^ instead of DFT+ [/, which usually produces 
adequate energetics for uranium oxides even though the predicted electronic structure might be wrong. ^ The validity 
of such treatment is guaranteed by comparing the calculated lattice parameters with experimental data. The cutoff for 
kinetic energy of plane waves was set as high as 500 eV to converge the electronic energy within a few rneV per atom. 
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Integrations over reciprocal space were performed in the irreducible Brillouin zone with about 8^-^36 non-equivalent 
k-points, depending on the system size. 

In this work, we did not take spin-orbit coupling into account. Previous calculations on metallic a-U,— perfect 
U02,^ and UsOg'*^ indicated that spin-orbit coupling has little effect on the structure and energy variation of uranium 
oxides, although it is well known that it is remarkable for magnetic properties. The only exception was reported by Yu 
et alM- Because inclusion of spin-orbit coupling in DFT-f- U deteriorates the already notorious problem of metastable 
states and convergence, and they did not employ any special technique to tackle this problem, it is highly possible 
that the reported results are due to premature calculations that trapped in metastable states. This can be inferred 
from the jagged variation (which is quite unphysical) of the calculated lattice constant with U, as shown in Fig. 3 of 
Ref. All of this indicates that their conclusion is still inconclusive. Based on these considerations, we neglected 
spin-orbit coupling in all of the following calculations. 

The supercell method has been employed with periodic boundary conditions. The arrangement of the structures is 
listed in Table HI in which Cm„ denotes a defective UO2 structure consisting of m fluorite cubic cells with n oxygen 
interstitials (or vacancies if n is negative)!^ The superscript u indicates the defect is with respect to uranium. For 
example, configuration "C81 has the same geometry as C81 but replacing the interstitial oxygen with a uranium atom, 
and "C8_i or C8_i corresponds to removing one lattice atom (uranium and oxygen, respectively) from a structure 
with 8 fluorite cubic cells (2x2x2). Before calculating the energy, all structures were fully relaxed until residual 
Hellman-Feynman forces (stress) were less than 0.01 eV/A. 

Cohesive energies at different volumes are extracted from the total energies by subtracting the contribution of 
isolated spin-polarized atoms. They are then fitted to a Birch-Murnaghan equation of state, 
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(3) 



to facilitate post-analysis, where the equilibrium atomic volume Vq can be determined from the effective lattice 
parameter as listed in Table HI With this equation, the pressure at K can be expressed as 



Pc{V) = -Bo 



^ "1 



(4) 



B. Formation Gibbs free energy 



Phonon contribution to free energy and thermodynamic properties can be computed via lattice dynamics or molec- 
ular dynamics simulations. After getting the vibrational density of states g{uj), the thermal part of the Helmholtz 
free energy per atom in quasi-harmonic approximation reads 

/•oo 

Fph{V,T)^kBT g{uj)\n[2smh{huj/2kBT)] doj. (5) 
"'0 

This method is accurate when anharmonic effect is insignificant. But the calculation is extremely demanding for 
ab initio simulations, especially for defective structures containing more than 90 atoms without any symmetry, as 
discussed here. Crude approximation includes computing the inter-atomic forces from empirical potentials'^ instead 
of by DFT calculations, while for off-stoichiometric compounds it is difficult to construct an accurate inter-atomic 
potential that maintains the charge neutrality of the cell generally. 

For practical purposes, we resorted to the Debye model, which is less accurate but much easier for computation. 
More important is that it provides a convenient description of lattice dynamics with a qualitative picture (even 
quantitative sometimes), and this is enough for our discussion here. In the Debye model, Eq.([5l) reduces to 



Fph{V, T) = SfcsTln [1 - e^^i^Qo/T)] - ksTf (Qd/T) + IksOD 



(6) 



where ks is the Boltzmann constant and / is the Debye function. The Debye temperature is approximated ai 
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1/3 



.1 1/2 



(7) 



Here the superscript p denotes the reference state (here the defect-free UO2), B is the bulk modulus, V is the effective 
volume per atom, and AI is the effective atomic weight. This approximation works well for a lot of materials^i^ 
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and is a physically sound extrapolation for 8d of U02±a; from that of perfect UO2. The value of the latter, i.e. Q^, 
takes 383 K at ambient conditions according to recent x-ray diffraction analysisi ^"'^^ 
The Gibbs free energy for a specific defective structure i is then given by 

G,(P, T) = E,iV) + FphiV, T) + PV. (8) 

For intrinsic point defects, the formation Gibbs free energy (FGE) of a Frenkel pair (FP) of species X is expressed as 

Gx.FP = Gl-^ + G^+^-2G^, (9) 

and for the Schottky defect (S) as 

Gs = G^;' + 2GZ;' - 3^G^- (10) 

If oxygen clusters are considered, we have to introduce isolated oxygen vacancy (Ou) to maintain the composition 
equation when in the closed regime where no particle exchange with the exterior occurs. For cuboctahedral clusters 
(COT), we thus should consider the Ot,-GOT complex and its FGE reads 

GcoT = G^+^ + yG^;i - [y + 1)G^, (11) 

where y is the excess number of oxygen and takes 5 in the COT-o clusterj^i^l In these equations, N is the number of 
atoms in a defect-free cell (96 for a 2 x 2 x 2 supercell) and is the corresponding Gibbs free energy. G^^^. and 

^COT ^^"^ ^^'^ Gibbs free energy of the cell containing the respective defect. 

It is worthwhile to mention that the COT cluster was originally proposed by experimentalists to explain the 
measured neutron diffraction patterns of \]40^^^ and XJ^Ot^^^ Historically, it was denoted by MgXag (or M6X37 
if an additional anion occupies the center) i^^i^^ However, it is a little misleading since the actual defect is composed 
of only 12 interstitials (forming the cuboctahedral geometry) and eight oxygen vacancies. There are two possible 
configurations for COT: with its center being occupied by another oxygen (COT-o) or being empty (COT-w). It was 
suggested that COT clusters should also present in \J02+x^^'^^^ Its geometry, energetics, and structural stability 
have been discussed in Ref.[27,|]. 

III. RESULTS AND DISCUSSION 
A. Properties at zero Kelvin 

1. Uranium monoxide 

In order to understand the general oxidation effects on compression behavior of uranium oxides, as well as to evaluate 
the significance of point defects on EOS of U02±a;, we studied UO, the a phase of U3O8, and the S phase of UO3. UO 
is face-centered cubic with the sodium chloride-type structure. The experimental lattice parameter is 4.92 ± 0.02 A.— 
However, there have been no claims of bulk UO having ever been synthesized, and reports of a UO surface phase 
on U metal for low exposures to oxygen^^ and UO thin films^^ could not be reproduced.^* Since UO has never been 
prepared with great purity, it has been suggested that the observed thin films actually represent uranium oxynitrides 
and oxycarbides, which form in the presence of N2 or C, and at low oxygen pressure<^ The equilibrium lattice 
parameter and electronic structure of UO have been studied theoretically by density functional method^^i^ but no 
mechanical property and compression behavior were explored, which are important for the purpose of understanding 
the response of the materials being oxidized from metal to dioxide and then to trioxide. 

Our calculated electronic density of states (DOS) is shown in the upper panel of FiglT] For UO, the /-electron 
lingers between the localized and delocalized pictures. Here we assumed itinerant / electrons and in the ground state 
it is metallic, as can be seen from the nonzero DOS at the Fermi level. The O atom has two unoccupied p states, 
and in the monoxide, the p band can accommodate two electrons from the uranium atoms (7s^) through charge 
transfer and hybridize with 6p^5/^6d^ weakly. The remaining valence electrons, including 6d and partially localized 
5/ electrons, start filling the conduction band, with the Fermi level pinned to the narrow / band. This agrees with 
earlier self-interaction-corrected local spin-density (SIC-LSD) approximation calculations.''^ The energy curves shown 
in Figl2] confirm that UO is in ferromagnetic (FM) ordering, with the anti-ferromagnetic (AFM) and non-magnetic 
(NM) configurations having an energy per cubic cell about 0.43 and 0.72 eV higher, respectively. A NM ordering state 
gives a shorter lattice parameter (4.84 A) and a higher bulk modulus (201 GPa), compared with early LDA results of 
4.88 A and 198 GPa^^ Meanwhile AFM and FM ordering have almost the same lattice parameter of 4.88 A, which is 
in good agreement with the experimental^^ value of 4.92 A and the SIC-LSD^^ result of 4.94 A. 
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FIG. 1. (Color online) The total and partial electronic density of states that projected onto the U5/ and O 2p orbitals for UO 
(upper panel), a-UsOg (middle panel), and (S-UO3 (lower panel) at OGPa. Vertical dotted lines indicate the position of the 
Fermi levels. 



TABLE II. Calculated lattice parameters (in A) and cell volume (in A^) of a-UaOg compared with other theoretical and 
experimental results. 



Approach 



Vol 



PBE 

PEE 

GGA+ 

Expt.^ 

Expt..^ 



FM, half-metaUic 
NM, metallic 
FM, insulator 
Insulator 
Insulator 



4.16 
4.16 
4.21 
4.15 
4.14 



11.57 
11.77 
11.61 
11.95 
11.96 



7.03 
6.87 
7.20 
6.72 
6.72 



338.36 
336.38 
351.92 
333.26 
332.74 



2. Triuramum octoxide and uranium trioxide 



UaOg is present as an olive green to black, odorless solid. In the presence of oxygen, uranium dioxide is oxidized to 
U3O8, whereas UO3 loses oxygen at temperatures above 773 K and is reduced to U3O8.— It is generally considered 
to be the more attractive form for disposal purposes of nuclear fuel because, under normal environmental conditions, 
U3O8 is one of the most kinetically and thermodynamically stable forms of uranium and also because it is the form of 
uranium found in nature. At low temperature U3O8 is in the a phase, a layered orthorhombic structure with space 
group Amr7i2,^-'*'- where the layers are bridged by oxygen atoms. Each layer contains uranium atoms which are in 
different coordination environments. 

We studied the a-U308 phase by computing its atomic geometry, electronic structure, and energetics. Table [TT] 
lists the calculated structure parameters by comparing with GGA-f C/^^ and experimental^^^ results. It can be seen 
that although pure GGA calculations with PBE functional fail to predict the correct insulator state, the optimized 
geometry is in good agreement with experimental measurements, except a slightly shorter b axis and a longer c axis 
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FIG. 2. Total energy variation of UO in different magnetic ordering states as a function of the cubic lattice parameter. 



were obtained. Deviation of the predicted equilibrium volume from the measured ones is less than 2%. In contrast, 
GGA-f U predicts a much worse result, with the equilibrium volume about 6% larger.*^ This is quite similar to the 
embarrassing situation encountered by the GGA+U approach in UO2, where a remarkable underestimation of the 
structural binding was also obtainedi^ Considering the acceptable performance of pure GGA on the energetics and 
EOS of UO2, we can expect that it will give reasonable results for the properties of U3O8 of interest here, even if a 
wrong half-metallic groundstate was predicted. 

The electronic structure of non- magnetic a-UsOg shows that it is metallic with the Fermi level pinned at the bottom 
of the conduction band. The energy gap below it is about 2.2 eV. The FM ordering phase has an energy about 0.3 eV 
lower per conventional cell than the NM ordering one, and is half-metallic. This energy gain is from pushing the 
minor spin component of the 5/ band away from the Fermi level. However, the remaining 5/ electrons on uranium 
atoms occupy the major spin band and partially localize at the Fermi level (see the middle panel of FigdJ . If Hubbard 
on-site correction were applied, the narrow 5/ band would be completely pushed away from the Fermi level and result 
in an insulator solution. 

Local magnetic moment is a helpful indicator for the localization degree of 5/ electrons. It is about 2 fiB per 
uranium atom for UO and 1.93 /is for U02,^ suggesting there are almost two 5/ electrons being partially localized. 
The difference is that in UO the localized electrons are at the Fermi level in the conduction band, whereas in UO2 
the on-site Hubbard interaction splits the 5/ band and leads to an energy gap (i.e., being Mott insulator), thus the 
electrons are localized in the valence band. On the other hand, UaOs is close to being a charge-transfer insulator, 
in which 5/ electrons are absorbed by oxygen to fill the empty 2p orbitals. But this process is incomplete and the 
remaining partially localized 5/ electrons in the conduction band pinned at the Fermi level contribute about 0.7 fiB 
per uranium atom, and make UsOg a half-metallic ferromagnet in GGA approximation. To open up the energy gap, 
the on-site Coulomb repulsion correction must be involved, and thus UaOg is an intermediate phase between the Mott 
insulator and the charge transfer insulator, as shown by the GGA-I- U calculatioujiS, This picture is in accordance 
with nominal chemical valence counting: on average each uranium atom bonds with 0.67 more oxygen in UaOg than 
in UO2, thus it loses about 1.3 more electrons in the former compound. Namely, most parts of the two localized 5/ 
electrons has been transferred to oxygen and only a small portion (0.7 electrons) is left behind. This charge transfer 
naturally creates a sharp energy gap just below the Fermi level, resulting in an almost closed shell on the uranium 
cores. We therefore predict that UO3 would be a pure charge-transfer insulator with the bare uranium ions in the 
closed-shell Rn configuration. 

Calculations on the 6 phase of UO3 confirms this argument. (J-UOa is a cubic phase in ReOa structure with a space 
group Pm3m^ The lower panel in FiglT] illustrates its calculated electronic DOS. It is a non-magnetic insulator. No 
any local magnetic moment on the uranium atoms can be observed. A distinct band gap of 1.6 eV can be clearly 
seen, which is due to the transfer of electrons to oxygen atoms, and leaves behind an empty atomic like 5/ band. As 
a result, the Fermi level moves to the top of the valence band. This completes the electronic structure evolution of 
uranium oxides along the oxidation process, i.e., from metal (a-U and UO, NM/FM) to Mott insulator (UO2, AFM), 
to intermediate insulator (U3O8, FM), and then to pure charge transfer insulator (UO3, NM) with a gradual shifting 
of the 5/ electrons from uranium to oxygen atoms (see FigH]). The correlation effect among 5/ orbitals also changes 
from one end to the other: in UO2 it is the strongest, whereas in UO3 we expect no strongly correlated effect to 
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FIG. 3. (Color online) Compression curves of UO, Q-UsOg, and iJ-UOs at zero Kelvin compared with that of U02±a;; the latter 
spans a quite broad pressure range. 



present at all. This can be verified by the correct prediction of the insulator groundstate of J-UOs with the pure GGA 
approach, as well as the excellent agreement of the cubic lattice parameter between the calculated value of 4.168 A 
and the experimental one of 4.165 A.^ 



3. Compressibility of U02±x 



The ordered phase of U4O9 has a complicated structure consisting of 4 x 4 x 4 fluorite-type subcells,'^'^ which is 
inaccessible for direct ab initio modeling because of the huge computational demands. We thus turn to the study of 
the Cli structure, as well as other defective U02±£c configurations as listed in Table HI and hope they can provide 
useful information on the general trend of oxidation effects on compression behavior. In Cli configuration there is one 
oxygen interstitial occupying the largest octahedral hole (i, i, i), and has the same chemical composition as U4O9. 
But one should keep in mind that it is more a defective UO2 rather than being U4O9, and would overestimate the 
bulk modulus slightly because of its high symmetry. 

The OK compression curves of defective U02±a: are shown in Fig|3l together with that of UO, a-UaOs, and J-UOa. 
The corresponding parameters for these curves are listed in Table HI It can be seen that although there might be 
some structure-dependence, a general trend is evident. With an increase of the oxidation degree from UO to UO2, 
the compound becomes harder and harder to compress, with the Cli configuration of U02±2: being the stiffest one. 
Further oxidization softens the material, and 5-UO3 is the most compressible one. Please note that a-U^Os and J-UOa 
transform into other structures under compression at a pressure scale of several GPa. This makes the P-V lines kink 
at the transition points, and increase the compressibility further. A detailed analysis of these pressure-induced phase 
transitions are, however, beyond the scope of this paper. 

The tendency due to oxidation also clearly manifests in FigHl where the correlation between cohesive energy, bulk 
modulus, and atomic volume is striking. Note that the values of a-U are taken from experiments,—'^ and others 
are the calculated results as listed in Table HI We can see that stoichiometric UO2 has the largest cohesive energy 
per atom, thus being the most energetically stable phase, and reflects the fact that when heated to 1573 K U3O8 
decomposes to UO2. As mentioned above, UO2, as well as U02±a;, lies in the intermediate stage of the charge transfer 
process. Therefore it has the strongest covalent bonds that arise from hybridization between U 5/ and O 2p orbitals, 
which then give rise to the smallest atomic volume, the deepest cohesive energy, and the highest bulk modulus, and 
therefore the least compressibility of U02±a; as shown in Fig|31 Further oxidization to U3O8 shifts the material quickly 
into ionic character, and the covalent bonds become weak. This leads to a huge volume expansion of more than 30%, 
and is dangerous for most nuclear application purposes. 

It is worthwhile to point out that in FigOthe P-V lines of U02±£c span a wide range of pressure, and is comparable 
to the difference between stoichiometric UO2 and other oxides. This underlines the importance of an accurate non- 
stoichiometric EOS, and it is inappropriate to represent it by any specific defective configuration, especially when the 
defect species change with temperature and pressure. Therefore a full EOS model is required for understanding the 
general thermophysical properties of off-stoichiometric compounds. That is, one has to study how the stoichiometry 
variation is accommodated in the material, what kind of microscopic entities are they in physics, and how these 
entities respond to different physical conditions. In the case of UO2, it is well known that point defects as well as 
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FIG. 4. (Color online) Variations of the cohesive energy per atom (D), the bulk modulus (Bo), and the atomic volume (Vb) as 
a function of composition deviation a; at GPa. 

some oxygen clusters account for the stoichiometry variation while keeping the fluorite-type matrix almost unchanged. 
These features make the construction of the non-stoichiometric EOS much easier. 



B. Finite pressures and temperatures 

Almost all of the earlier defect population analyses were carried out with K formation energies i^i^i^ i^^i^^'^^ It is 
widely assumed that the contribution from lattice vibrations is insignificant, but lacks detailed assessment. Also local 
residual stress and thermal expansion effects were completely neglected. The following discussion is devoted to this 
issue, and reveals that thermal vibrations have a limited but important influence on defect populations, whereas the 
pressure effect is much more striking in comparison. 

It is worthwhile to note that since UO2 is a semiconductor, the defects can be charged i^"'^^ These charged species 
might have several charge states^i^ These facts and the slow convergence of a charged system in computation com- 
plicate the discussion (the energy converges linearly with respect to the inverse of the supercell size at the first order 
correction). Since our main purpose here is to investigate the phonon contribution and the compression effects on de- 
fects behavior, which should be similar for neutral and charged defects (we are reminded of the qualitative consistency 
of these two kinds of defects on OK formation energyiSiii^) , we thus focus our discussion here on neutral defects only. 
Extension of the analysis to charged defects is straightforward as long as the energetics of these species are known. 



1. Variation of formation Gibbs free energy (FGE) 

Variations of FGE as a function of temperature and pressure are shown in FigjS] for intrinsic point defects (oxygen 
and uranium Frenkel pairs and Schottky defect, respectively) and in FiglHlfor the 0„-COT-o complex, as defined in 
Ea.([TT|). It can be seen that phonon contribution reduces the FGE of the uranium Frenkel pair (FP), whereas it 
enhances that of the oxygen FP. As a result, its contribution to Schottky defect almost completely disappears due 
to cancellation of the opposite effect on uranium and oxygen defects. The maximum difference of the FGE between 
and 1800 K is less than 0.62 eV for Frenkel pairs and 0.12 eV for Schottky defect, respectively. In other words, it 
counts only about 11% for oxygen FP, 4% for uranium FP, and 1% for Schottky defect. 

The calculated FGE at GPa and K are consistent with that of Ref. Q . But we have to point out that there is 
a typo in Ref.^] where the formation energy of the uraniumpoint interstitial (U-Int.) should be 5.642 eV instead 
of 8.194 eV. This change affects Table I and Table IV in Ref.Q. Also, it reduces the energy of the uranium Frenkel 
pair (U-FP) in Table IV from 17.2 to 14.6 eV. However, this does not affect the defect population analysis since the 
concentration of uranium interstitial is too small to have any significant contribution to the stoichiometry of UO2 . 

On the other hand, lattice vibrations decrease the FGE of the COT-o complex (see FiglH]). Increasing the tem- 
perature from to 1800 K, the change in the FGE of the COT-o complex is less than 0.72 eV, which is about 3% of 
the total formation energy. Compression depresses the phonon effects slightly. But totally, the thermal vibrational 
contribution to the FGE of point defects and the oxygen defect cluster is small. 

In contrast, the compression effect is quite prominent. Within a pressure range of 30 GPa, which is easy to achieve 
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FIG. 5. (Color online) Variation of the formation Gibbs free energy of Frenkel pairs and Schottky defect in UO2 under 
compression at different temperatures. Arrows indicate the direction of temperature increase with a step size of 600 K. Dashed 
lines correspond to K results with zero point energy included. 
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FIG. 6. (Color online) Formation Gibbs free energy of the COT-o complex as a function of temperature and pressure. 



with diamond-anvil cell (DAC) devices or micro shock-wave loadings in nuclear materials driven by thermal spikes 
or collision cascade events arisen from fission reactions, the FGE increases 1.07, —1.23, and 4.02 eV for uranium FP, 
oxygen FP, and Schottky defects, and amounts to 7%, —23%, and 38% of the zero pressure values, respectively. For 
the COT-o complex, the FGE change is — 4.37eV, about —17% of the OGPa value. Note that a pressure of 30 GPa 
corresponds to an energy scale of 2.4 eV per atom in UO2. 



2. Defect populations 

In addition to being a key content in the construction of a non-stoichiometric EOSj^ the importance of analyzing 
defect populations, or equivalently the concentrations, also lies in the challenge of the elusive physics of U02+a:. 
Although experiments have shown that COT-o clusters constitute the ordered phase of U4O9 on a fluorite matrix,—"— 
it is still in debate whether or not such a cluster also appears in U02+x'^'^'^'^'^ Calculations using OK formation 
energy suggested that the existence of COT-o in UO2+2; is quite likely, but the observed ratio of oxygen interstitials 
occupying different sites cannot be accounted for quantitatively.^"'-™ This implies that we cannot exclude other 
possibilities. In this sense it is very important to figure out how the population of each possible defect changes 
with temperature and pressure. On the other hand, as mentioned above, one of the GGA-I- U calculations predicted 
a positive OK formation energy for point oxygen interstitial.— If this reflects the true physical behavior, we can 
expect that phonon contribution at elevated temperatures should be as such that it makes oxygen defects the major 
component. These considerations underline the necessity to investigate the defect populations at different pressures 
and temperatures with phonon contributions having been taken into account. 
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FIG. 7. Variation of concentrations of point defects at a temperature of 1700 K with x, the deviation from the stoichiometry 
of UO2 in the hyper-stoichiometric regime: left is at a pressure of GPa and right is under 10 GPa. Solid (dotted and dashed) 
line indicates the concentration in oxygen interstitial (oxygen vacancy and uranium vacancy). The concentration of uranium 
interstitial is negligible. 
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FIG. 8. The influence of oxygen cluster COT-o on defect concentrations at and 10 GPa, respectively. The dash-dotted line 
indicates the concentration of COT-o cluster; others are the same as in FiglT] 



Using the point defect model (PDM) developed by Lidiard and Matzke^ i^^i^^'^^ defect populations around the 
dilute limit can be computed. The key point is that here the FGE instead of the K formation energies are employed. 
In the following we mainly discuss the hyperstoichiometric {x > 0) regime because in the other side of a; < 0, it is 
always oxygen vacancy (Oy) that predominates. Although both pressure and temperature enhance the concentration 
of point oxygen interstitial (O^) and uranium interstitial (Ui) slightly, they are still at least two orders smaller than 
that of Ov and can be neglected safely. 

Figure [7] illustrates the defect populations as a function of composition at 1700 K under pressures of and 10 GPa, 
respectively. By comparing with Fig. 11 in Ref. '9^ where only K and GPa formation energies were used in the PDM 
analysis, we find a slight decrease of the concentration of 0^ and an increase of the uranium vacancy (U„), and the 
latter becomes the major component when x > 10^'^. Both lattice vibrations and thermal expansion contribute to 
this change. From FigOwe learn that phonon favors uranium defects. Thermal expansion also prefers uranium FP 
and Schottky defect because at an elevated temperature thermal expansion leads to a larger volume that corresponds 
to a negative cold pressure, which reduces the formation energy of uranium FP and Schottky but lifts that of oxygen 
FP. This result is strongly against the experimental observation that oxygen defects are predominant in UO2 at any 
temperature and stoichiometry i^'^i^^'^^ Increasing pressure can depress the uranium vacancy, as shown in the right 
panel of FiglT) This is easy to understand by inspecting the variation of FGE with pressure as shown in FigIS] 
Nevertheless, this compression effect is not enough to explain the observed predominance of oxygen defects because 
the possible residual stress within the fuel is usually less than 1 GPa, and at such a pressure scale the uranium vacancy 
is still comparable to oxygen defects. 

It thus becomes interesting to investigate whether phonon and compression effects can alter the behavior of the 
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FIG. 9. The same as in FigObut now with a temperature of 700 K. 



COT-o cluster greatly or not. Earlier studies showed that when only OK formation energies were used, the COT- 
o cluster is dominant at low temperature and high composition region.-^^'^^ Figure |8] plots the variation of defect 
populations at 1700 K with COT-o included. Note here that we didnot consider the internal degeneracy of this 
cluster," since physically no two clusters can occupy the same lattice site simultaneously if we ignore quantum effects 
and treat oxygen as a classical particle. We can see that, as in the previous case, phonon contribution enhances 
uranium vacancy, which is very clear if comparing with Fig. 3 of Ref.[27]. The earlier conclusion that COT-o is 
predominant is still preserved when x > 10~^. At a pressure of lOGPa, uranium vacancy is drastically suppressed, 
leaving only single oxygen interstitial and COT-o in competition (the right panel of FiglS]). Lowering the temperature 
to 700 K further suppresses uranium vacancy and oxygen interstitial, as Figl9] shows, meanwhile COT-o gets some 
promotion. Increasing the pressure to 10 GPa then sweeps all other defects away and leaves only Oi and COT-o 
clusters, with a slight depression of the latter and an enhancement of the former. 

It is evident that lattice vibrations and hydrostatic compression change the competition among defects greatly, and 
the conclusions obtained with K and GPa formation energies need be revised: when considering only point defects, 
the earlier predicted predominance of 0^ doesnot hold anymore with the full FGE, and \Jy becomes the dominant 
defect when x > 10~^. That is to say, at the full level of the theory the experimental fact cannot be reproduced if 
only neutral point defects are considered. On the other hand, when the COT-o cluster is taken into account, the 
predominance of this oxygen cluster can be easily established both at the full FGE level and at the K and GPa 
formation energy level. The calculated results clearly demonstrate that hydrostatic pressure has a strong influence 
on defect populations, and oxygen clusters are definitely necessary in order to suppress the concentration of uranium 
vacancy to match the experimental observation qualitatively. 

Considering the low temperature phase diagram^^ of U02+a; and the negative formation free energy of 0^ and COT- 
o, it becomes clear that the possibility of phase separation between UO2 and U4O9 must be considered seriously if one 
wants to understand the underlying physics correctly. Because U4O9 can be viewed as an ordered phase of COT-o 
clusters on a fiuorite matrixj^i^^ this phase separation is in fact a process of formation and segregation of oxygen 
clusters. Thus at the intermediate composition and elevated temperature range where U02+a; is stable, the dynamics 
of decomposition/formation of COT-o is crucial, as well as other competitive clusters such as {V-30 )2fi^i^'2^i^ It 
is worthwhile to highlight the connection between V-30 based clusterai^'^'^ and cuboctahedral cluster o^^'^^'^^i^^ 
(mainly COT-o). Although the latter is more stable at low temperature it was predicted that with increasing 
temperature COT-o will gradually dissociate into a set of single point defects<^ The intermediate state of this process 
is of great importance for understanding the kinetic behavior. The relatively low energy of (1^-30 )„ makes it a 
promising candidate of this transition state^Si This is possible because the COT cluster in fact can be viewed as a 
complex formed by eight tightly connected 1^-30 clusters (a distorted version of V-30 ). Then partial dissociation 
of COT clusters will likely go through a configuration of V-30 /O kinetically to lower the energy barrier; that is. 



COT-o ^ {V-30")n ^ 50,. 
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3. N on- stoichiometric equation of state and shock Hugoniot 



OfF-stoichiometry is widely observed in materials of many categories, for example, in ionic compounds {e.g., oxides, 
nitrides, and carbides), alloys and intermetallic compounds, doped semiconductors, and high-temperature supercon- 
ductors, and so on. In these materials, the underlying structure of the matrix is usually unchanged or just slightly 
distorted, but the electronic and magnetic properties can be controlled and tuned for application purposes. The 
maximal stoichiometry deviation indicates the capability of the matrix structure to tolerate distortions caused by 
doped elements or defects. The mechanical properties usually are determined by the matrix, but are also varied with 
respect to different chemical compositions. 

The thermodynamic behavior of these materials across a wide range of temperature and pressure are of interest 
not only for fundamental physics and material science but also for engineering applications, in which EOS is one 
of the most important contents. There is an elegant theory developed on Ising models of lattice gases, namely, the 
cluster variation method (CVM)f2S which describes the phase diagram and EOS of binary and pseudo-binary alloys 
very wel L^^'^^" — To generalize this method beyond the ternary system is difficult. For ionic compounds, however, 
the stoichiometry usually is accommodated via defects, and is a multi-component system with the number of species 
greater than four. Furthermore, not all of these defects sit on the lattice site of the matrix, which complicates the 
construction of the theory, and slows down the convergence even if a CVM-like theory could be constructed.™ From this 
point of view, a different route should be attempted for a practical ab initio EOS model (i.e., without any empirical 
parameters) of non-stoichiometric compounds. 

Having known the defect concentrations, we can express the total Gibbs free energy G of a non-stoichiometric 
material at pressure P and temperature T as a function of the defect concentrations that accommodate the 
stoichiometry deviation. Around the vicinity of the dilute limit, Taylor expansion of G with respect to rii reads^- 



2 ^ 

where G^ is the Gibbs free energy of the matrix, Ai is the contribution of single defects, and Btj is the interaction 
strength between defects, and so on. At the first order approximation we can express the non-stoichiometric Gibbs 
free energy by making use of specific defective structures as reference; that is, 

G-G^ + E^-^- (13) 

Here the superscript "ref" indicates the defective reference system, rii is the concentration of defect type i, and the 
Gibbs free energy difference is AG^ = Gf^ — G^. Obviously this linear approximation assumes that the interaction 
term between defects 5^ is insignificant, which is only true at the dilute limit. However, making use of reference 
states extends the applicability of this approximation greatly. It is evident from Eq. ([13|) that G is close to being exact 
when Ui = nf^ . Nevertheless, within this linear approximation, the solubility of defects and spinodal decomposition 
of the system cannot be correctly described, where the nonlinear terms in Eq. (1121) are essential. 
Within this framework, the system volume becomes 

V = Vr> + J2 ' ,ef (14) 

and this gives the equation of state for a non-stoichiometric material. Defect concentrations at thermodynamic 
equilibrium conditions are determined with the composition equation using the point defect model as discussed above. 

Hugoniot is a set of locus of thermodynamic states that describe the response of a material to shock wave loadings. 
It is governed by the law of energy conservation across the discontinuous zone generated by shock waves, i.e., 

E-E^ = \^{P + P^) (Vb - V) , (15) 

and serves as an important content of high-pressure high-temperature EOS. For theoretical analysis, it is convenient 
to reformulate Ea. (fT5)l as^S 



G - Go = ^ (P - Po) 



dG\ fdGa 



where Go is the system Gibbs free energy at an initial condition of pressure Po E^nd temperature To. 
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FIG. 10. Hugoniot of U02+i: shocked from an initial state of 500 K and QGPa. 

Taking the largest defective configurations as the reference structures {i.e., C8 and COT-o configurations as fisted 
in Table U that are based on a 2 x 2 x 2 supercell of fluorite cubes), the shock Hugoniot of \]02+x up to 70 GPa were 
calculated and shown in FiglTU]with x = 0, 0.1, and 0.25, respectively. Note that since here all defects are defined 
on a fluorite cubic cell of UO2, the Gibbs free energy G must be given with respect to the same cell. The initial 
condition of shock is at OGPa and 500 K. The highest shock temperature is less than 1200 K, reflecting the high bulk 
modulus of the material and thus a weak energy accumulation of shock waves. Although the cubic fluorite phase of 
UO2 transforms into orthorhombic Pnma at room temperature beyond 42 GPa,^'^ the transition is rather sluggish 
and the fluorite phase persists up to 70 GPa. ^ This suggests the selected range of pressure in FigHniis appropriate. 
Beyond that pressure, we must reanalyze the geometry of defects and their occupation dynamics because the matrix 
has already been changed due to phase transition. 

From Figlin] we can see that the Hugoniot of 1102+3; is very similar to stoichiometric UO2 , with the whole curve 
shifting to the left and reducing the cell volume correspondingly with an increase of x. This is reasonable because 
from the above analysis we learned that within this pressure and temperature range the COT-o cluster is always 
predominant, and the defective structure with this cluster has a very similar OK energy curve as that of UO2 (see 
Table H]). However, this is not always the case for non-stoichiometric compounds. Along the shock wave loading path, 
pressure and temperature change from point to point because of energy accumulations. If species of the major defect 
component were changed due to this kind of thermodynamic condition variation, the atomic scale behavior of defects 
would then manifest itself on the Hugoniot curve. In that case the non-stoichiometric EOS would be a useful tool to 
study this phenomenon theoretically and thoroughly. 

IV. CONCLUSION 

We have presented calculations, using the LSDA-1- U exchange-correlation functional, of the high-temperature high- 
pressure properties of point defects in U02±£!;. The phonon contributions were taken into account via the Debye 
model. The results showed that lattice vibrations enhance the population of uranium vacancy slightly, whereas 
hydrostatic compression changes the energetics of all defects and then their populations greatly. At ambient pressure, 
however, inclusion of phonon free energy and thermal expansion effect fails to reproduce the experimental observation 
that oxygen defects should be predominant at the x > regime, suggesting that the inclusion of defect clusters 
such as COT-o is necessary. Furthermore, both lattice vibrations and compression favor the COT-o cluster, and 
its predominance in the hyper-stoichiometric regime gets reinforced as compared with the results obtained with K 
formation energies. 

In order to understand the overall trend of oxidation effects, the K electronic structures and equations of state of 
UO, U02±a:, a-UsOg, and (J-UOa were investigated thoroughly, revealing the transition of the groundstate of uranium 
oxides from being metallic to Mott insulator and then to charge-transfer insulator across the full oxidation degree. 
Variations of cohesive energy, atomic volume, and stiffness of the oxides due to evolution of the chemical bonds from 
being metallic to covalent and then to ionic characteristics were well demonstrated. The necessity to develop an EOS 
model for off-stoichiometric compounds was also highlighted by comparing the compression curves of defective U02±a: 
with that of other uranium oxides. With the knowledge of defect populations within a wide range of pressure and 
temperature, a specific equation of state for non-stoichiometric compounds was constructed. The shock Hugoniots of 
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UO2+X were then studied to illustrate the performance of this model. 
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